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CORRECTIBLE ERROR PATTERNS 
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Abstract. We demonstrate a majority-logic decoding algorithm for de- 
jrt ' coding the generalised hyperoctahedral group Cm wr Sn when thought 

*^ of as an error-correcting code. We also find the complexity of this de- 

^—i coding algorithm and compare it with that of another, more general, 

|~ — algorithm. Finally, we enumerate the number of error patterns exceed- 

rvq ing the correction capability that can be successfully decoded by this 

algorithm, and analyse this asymptotically. 

\^ Keywords: Error-correcting code, permutation code, asymptotic enu- 

r) meration. 

^ MSC2010 classification: 94B35 (primary), 05A16, 20B05, 94B25 (sec- 
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1. Introduction 



,^j. Sets, or groups, of permutations may be used as error-correcting codes, 

^^ with permutations in list form as the codewords, and the usual Hamming 

[^ distance. This idea goes back to the 1970s, for instance to the papers of 

^^ Blake [3] and Blake, Cohen and Deza [3]. More recently, there has been 

OO a resurgence of interest in such codes because of a potential application 

T\ to "powerline communications", where electrical power cables are used to 

J>" transmit data as well as electricity. For instance, the 2004 paper by Chu, 

K^ Colbourn and Dukes [7] gives a description of this, and some constructions 

5-H for suitable codes, while a more general survey can be found in Huczynska's 

"^ 2006 paper [H]. 

In [2], the first author gives a decoding algorithm which works for arbi- 
trary permutation groups when used as codes in this way. In [Ij, he considers 
certain families of permutation groups in more detail. In order to describe 
the case we are interested in in this paper, we need the following background 
material. 

Suppose H and K are permutation groups acting on sets F and A respec- 
tively, where |F| = TTi and A = {1, . . . , n}. The wreath product G = Hwr K 
is constructed as follows. We consider the action of the Cartesian product 
H^ = HxHx---xH on n disjoint copies of the set F, labelled by the 
elements of A. We then form the semidirect product of iJ" with K, where 
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K acts on H"' according to its action on A = {1, . . . , n}; the resulting group 
is G = H WT K. Now, we can define an equivalence relation on F x A where 
the equivalence classes are the copies of F; this equivalence relation is pre- 
served by the action of G, and so forms a system of imprimitivity or block 
system for G. (See Cameron [B] for more information about permutation 
groups.) 

One family considered in |lj were the groups H wi Sn, where H isa regular 
permutation group of order m. In this paper, we consider the special case 
of this where H is a cyclic group of order m, so we have G = Cm wr Sn 
acting in its imprimitive action on n copies of {1, . . . ,m}. We call these 
groups generalised hyperoctahedral groups, as in the case m = 2 we have 
the well-known hyperoctahedral group (the automorphism group of the n- 
dimensional hypercube) . 

In Section [2] we give an alternative decoding algorithm from that given 
in [2j, that can only be used in this case, and in Section ^ we show that 
this algorithm is better-performing (in terms of time and space complexity) . 
Finally, in Section [4] we show that certain patterns of more than r errors 
can be successfully decoded by this algorithm, and in Section [5] we analyse 
the asymptotic behaviour of this. 

Recall that the minimum distance, d, of a code is the least value of the 
Hamming distance over all possible pairs of codewords, and that the cor- 
rection capability (i. e. the number of errors that can be guaranteed to be 
corrected), r, is given by r = [^^J. 

Proposition 1. The correction capability of G = Cm'^'c Sn is r = [^^^^^J. 

Proof. The Hamming distance between two permutations g and h is precisely 
n— \Fix{gh~^)\, where Fix(g') is the set of points fixed by that permutation. 
Using the group structure, the minimum distance of a permutation group G 
of degree n is therefore equal to 

n — max |Fix(5()| 
geG 

gf^id 

(see [2] for further details). Now suppose G = CmWY Sn- Each copy of 
{!,..., m} forms an imprimitivity block for G, and furthermore if an element 
of G fixes one point in a block, it must fix all m points in that block. 
Consequently, the maximum number of fixed points is {n — \)m (it is easy 
to construct elements with this many), so the minimum distance is nm — 
(n — l)m = m, and so the correction capability is r = [^^^^^J . □ 

2. The decoding algorithm 

The decoding algorithm we give below makes use of the relatively straight- 
forward combinatorial structure of the group G = Cm wr Sn ■ Since the 
group permutes n imprimitivity blocks of size m, a permutation in G when 
written in list form can be divided into n blocks of length m. The ordering 
of the blocks gives the action of Sn on the imprimitivity blocks, and the 
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(cyclic) ordering of the symbols within a block gives the cyclic permutation 
gives the corresponding element of Cm- Thus each position effectively holds 
two pieces of information which are constant throughout that block: a block 
label and a cyclic shift. 

Example 2. The following permutation is an element of C5 wr S4: 

[7, 8, 9, 10, 6 I 15, 11, 12, 13, 14 | 20, 16, 17, 18, 19 | 5, 1, 2, 3, 4]. 

As can be seen, the permutation splits into four blocks of length five, con- 
taining 1, . . . , 5, 6, . . . , 10, etc. 

Recall from Proposition 1 that the correction capability is r = [^^^^J- 
Consequently, if we assume there to be a maximum of r errors, there will be 
a majority of positions in each block which contain the correct symbol. The 
decoding algorithm uses this fact: the majority of the block labels and the 
majority of the cyclic shifts will be correct, so this allows the reconstruction 
of the transmitted word. 

If a list L has a unique most frequently-occurring element, we denote it 
by majority (L). 

Algorithm 3. Input the received word w = [wi\w2\ ■ ■ ■ \wn], where Wi = 
[wii, . . . , Wim]- For each i and j, we calculate qij and Sij where Wij = mqij + 
Sij, then map Wij to a pair (6jj, Cij) as follows: 



hij :-- 




if s,, + 
1 if Sij = 


Cij :-- 


= Sij-j 


mod m. 



Defining bi = majority[6ii, . . . , bim] and q = majority[cji, . . • , Qm] for each 
i, the list [61, ... , hn] gives a permutation of {0, 1, . . . , n — 1} corresponding 
to the element of Sn acting on the blocks, and the list [ci,...,cJi] gives 
the cyclic shifts within each block. We can then reconstruct the original 
permutation g := [gi\g2\ ■ ■ ■ \gn], where gi = [gn, . . .,gim], 9ij = mbi + Uj 
and tij = j + Ci mod m. 

Example 4. Suppose we transmit the following element g S C5 wr 84^: 
[7, 8, 9, 10, 6 I 15, 11, 12, 13, 14 | 20, 16, 17, 18, 19 | 5, 1, 2, 3, 4]. 

Then suppose we receive the following word w. 

[17,1,9,10,6 I 15,11,12,13,14 | 20,16,17,18,19 | 5,1,2,3,4]. 

This clearly has errors in positions 1 and 2. Having split this into four blocks 
of length five, we obtain the data shown in Table [!} 

Taking the "majority" elements, we find the block permutation /3 = 
[1, 2, 3, 0] and cyclic shifts of [1,4, 4, 4]. We have the information needed to 
reconstruct the transmitted permutation: for instance with i = and j = 1, 
we have 60 = 1; cq = 1) *oi = 1 + 1 mod 5 = 2 and so goi = 1x5 + 2 = 7. 
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i J 


Wij 


ibij,Cij) 




i 3 


Wij 


(&ii,%) 


1 


17 


(3,1) 


2 1 


20 


(3,4) 


2 


1 


(0,4) 




2 2 


16 


(3,4) 


3 


9 


(1,1) 




2 3 


17 


(3,4) 


4 


10 


(1,1) 




2 4 


18 


(3,4) 


5 


6 


(1,1) 




2 5 


19 


(3,4) 


1 1 


15 


(2,4) 


3 1 


5 


(0,4) 


1 2 


11 


(2,4) 




3 2 


1 


(0,4) 


1 3 


12 


(2,4) 




3 3 


2 


(0,4) 


1 4 


13 


(2,4) 




3 4 


3 


(0,4) 


1 5 


14 


(2,4) 




3 5 


4 


(0,4) 



Table 1. Data obtained during decoding in Example [4] 



Performing these calculations for each i and j, we can recover the transmit- 
ted permutation: 

[7, 8, 9, 10, 6 I 15, 11, 12, 13, 14 | 20, 16, 17, 18, 19 | 5, 1, 2, 3, 4]. 

We conclude this section by mentioning that the first author has imple- 
mented this algorithm in the computer algebra system GAP [9]. 



3. Complexity 

We recall that Algorithm [3] has three parts: calculating the numbers 
(bij,Cij), which involves integer arithmetic; determining the most frequently- 
occurring elements bi and q; and reconstructing the decoding permutation g, 
which involves more integer arithmetic. In order to determine the complexity 
of this algorithm, there are some assumptions we need to make first. 

• integer arithmetic can be done via a look-up table, in constant time; 

• comparing the sizes of two integers can be done in constant time; 

• finding position z in a list of length k takes 0(log A:) time. 

However, the second step is more complicated and requires the following 
lemma. 

Lemma 5. Let L be a list of length m with symbols chosen from S = 
{1, . . . , /c}. Suppose L has a unique most frequently occurring element, x G 
S. Then the time taken to determine x is 0{k + mlogk). 

Proof. We begin by producing an auxiliary list K of length k, initially set 
to [0, 0, . . . , 0]. This takes k units of time. We then work through each of 
the positions of L: in each position, we do as follows: 

• read the symbol, i, (taking one unit of time); 

• find position i in K (taking 0(logA;) time); 

• increment that entry by 1 (taking one unit of time). 
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This turns K into a list of the frequencies of each symbol in S in the list 
L. Doing this for each of the m entries of L requires a total of 0(?7ilogA;) 
time units. We then work through K to find the position of the maximum 
element. This will require 0{k) comparisons. Combining this, we have 
0(/c) + 0(m log k) + 0{k) = 0{k + m log k) as required. D 

We observe that the method described above is not necessarily the best 
possible; other methods may be faster, and which method is the best may 
depend on factors such as the relative sizes of m and k. However, we now 
use it to determine the time complexity of Algorithm [3j 

Theorem 6. The time required to perform Algorithm^is 0{mnlogm) (if 
m > n) or 0(n^ + rn-n log n) (if m < n). 

Proof. The first stage is the calculation of the numbers (6jj, Cij). There are 
mn such calculations to perform, and we have assumed that each takes a 
constant amount of time, requiring a total of 0{mn) time units. 

The next stage is, in each block i, to determine the most frequently occur- 
ring block label hi and most frequently occurring cyclic shift Ci. This involves 
determining the most frequently occurring element in a list of length m with 
symbols chosen from a set of size n (for the block labels) and from a list of 
length m with symbols chosen from a set of size m (for the cyclic shifts). 
By Lemma [5] above, the first of these will take 0(n + mlogn) time units, 
the second 0{m + mlogm). As this has to be done in each of n blocks, this 
gives a total of 0(n? + mn log n + mn + mn log m). 

The final stage is the reconstruction of g, which requires m integer arith- 
metic operations in each of the n blocks. As we have assumed that integer 
arithmetic takes constant time, this requires a total of 0{mn) time. So the 
total time required is O (mn) -|- O (n^ -|- mn log n -|- mn + mn log m) + O (mn) . 
If m, > n, this reduces to O (mn -|- mn log m) = O (mn log mn) , while if m, < n 
it reduces to 0(n^ -|- mnlogn), as required. D 

We should also consider the space complexity of Algorithm [3j This time, 
we require a look-up table for our integer arithmetic, and there are also 
items that have to be stored whilst the algorithm is being performed. 

Proposition 7. The amount of storage space required by the decoding algo- 
rithm is 0{mn'^), and the space required to perform the algorithm is 0{mn). 

Proof. We need to store a look-up table, where for each of the mn symbols, 
for n possible divisors we record a quotient/remainder pair. This requires 
a total of 2mn^ = 0{mn'^) storage units. To perform the algorithm, we 
need to store mn quotient /remainder pairs, then the two auxiliary lists (one 
of length m and one of length n) to find their most frequently-occurring 
element, and need mn units to store the reconstructed group element. This 
gives a total of 2mn + m-\- n + mn = 0{mn) units. D 

In [2j, a more general decoding algorithm was given, which works for 
arbitrary permutation groups; also in [2j its complexity was analysed in a 
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similar fashion to the above. In the case where the group is the generahsed 
hyperoctahedral group Cm wr Sn, bounds on the time and space complexity 
for the more general algorithm are bounded by 0{m?Ti?) and 0{'rrfirfi) re- 
spectively. Thus, from the point of view of a worst-case analysis, Theorem 
[G] and Proposition IT] suggest that Algorithm Is] is an improvement. (Both 
algorithms require 0{mn) space to perform the algorithm.) 

4. Enumerating correctible error patterns 

We observe that the algorithm described above only requires a majority 
of correct elements in each block. Consequently, there will be patterns of up 
to nr errors that can successfully be corrected, as long as there are no more 
than r errors in each block. In this section, we investigate how many such 
patterns there are. 

For positive integers k, n and r, define Vn,T{k) to be the set of all partitions 
of the integer k into at most n parts, and where each part has size at most 
r. For vr G Vn,r{k), we denote the number of parts of size i by /i(vr) (so that 
X] /i(7r) < n). We also define a quantity Ci(7r) to be 

i=i 

for i > 2, with ci(7r) = 0. That is, Ci(7r) is the number of parts in vr of size 
strictly less than i. 

Proposition 8. For a word in CmWi^ Sn, and for k < nr, the following 
number of patterns of k errors can be corrected: 



v^ j-r 1 1^ ~ CiiTT ) \ f m 



En,mAk)= Y. IIv f-(n) 



M^) 



Proof. For a pattern of k errors to be corrected, they can be spread across 
up to n blocks, as long as there are no more than r errors in each block. So 
for a given partition vr G Vn,r{,k), for each part of n we have to choose (i) 
which block contains that many errors and (ii) where in that block they lie. 
Working through i in increasing order, for each i there are n — Cj(7r) blocks 
remaining, of which we choose /i('/r) (corresponding to the /i(vr) parts of 
size i). Then in each of the /i(7r) blocks we have chosen, we choose i error 
positions from the m available. D 

Note that if /c > nr, the set Vn,r{k) is empty, so En,m,r{k) = 0. 

While this is a tidy combinatorial expression for the desired quantity 
En,m,r{k), its behaviour cannot easily be seen, especially as we wish to com- 
pare it with the total number of patterns of k errors, ("^"'). A first step 
would be to find a recurrence relation. 
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Lemma 9. The numbers En,m,r{k) satisfy the recurrence relation 

-l,m,r (,"' '')■ 



i^n,m,r\'i) — / ^ I , j J^n 
1=0 ^ ^ 



Proof. Suppose there are I errors in the n block; there are (™) ways of 
arranging these. Then there are k — I errors in the remaining n — 1 blocks, 
so there are En-i,m,rik — I) ways of arranging these. Summing over all 
possible values of / < r, we obtain the required relation. D 

This recurrence relation assists us in studying the generating function for 
En,m,r{k)- Let £n,m,r{x) denote this function, that is 

A;>0 

and define 

Proposition 10. The generating function £n,m,r{x) can be rewritten as 
Proof. Applying Lemma [9] and re-summing, we obtain 

fe>0 

r 



fc>0 1=0 ^ ^ 



By iterating this, and with the observation that £o.m,rix) = 1, we have 



\l=0 ^ ' ^ / 



as required. D 

The probability that a word with exactly k errors is correctible is then 
given by 

as the number of all possible patterns of k errors is given by ("™) , of which 
En,m,r{k) are correctible. (Note that the number of words related to a fixed 
pattern of k errors is {mn — 1)^, and that this factor is cancelled.) 

From the point of view of applications it is perhaps more useful to con- 
sider the probability Pn,m,r{p) that a word is correctible under the assump- 
tion that an individual error occurs with probability p. In other words, we 
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consider a probabilistic model in which the number k of errors is binomially 
distributed, k ~ B{mn,p), so that the expected number of errors is given 
by pmn. 

Proposition 11. For a word in Cm^^ Sn, if an individual error occurs 
with probability p then the probability that a pattern can be corrected is given 
by 

Pn,mAP) = (1 -P)""'?n,m,r {?/{!- p)) ■ 

Proof. A pattern is correctible if there are at most r errors in each block. 
The probability that exactly / errors occur in a block of size m is given by 
('J^)p'(l — p)"^~', so that the probability that there are at most r errors in 
each of n blocks is given by 



Pn,mAP)=(J2["l)p'i^-Py 



\m—l 
I / IP^'-P 
=0 

which gives the desired expression. D 

Of course these two probabilistic models are not equivalent. 

5. Asymptotic analysis 

5.1. Asymptotics of pn^m,r{k) for CmWi^ 5n as n — )• cxd. Let us focus on 
the asymptotics of pn,m,rik) when the error frequency k/mn is fixed. Here, 
we will discuss the case where n — t- oo as m is fixed. 

We first give an expression for -E'n,m.,r(^) which is amenable to asymptotic 
treatment. 

Lemma 12. 

En,mAk) = ^ / iFmAz)T ^fc^ > 

where the contour of integration is a counterclockwise circle about the origin. 



Proof. This follows directly from Proposition 10 and the Cauchy Integral 



Formula. D 



An asymptotic analysis of the integral in Lemma 12 is obtained from a 
saddle-point approximation. (See Flajolet and Sedgewick [8] for background 
material on this technique.) For n and k large, the behaviour of the integral 
is determined by the exponential of 

nlogFmA^) -klogz . 
There is a unique positive saddle Q given by 



-I 



k 

logFmAO logC 

n 



and the asymptotics are obtained by approximating the integrand around 
the saddle by a Gaussian. This is the content of Theorem VIII. 8 of [8j, 
which we use to obtain the following result. 
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Proposition 13. Let A = k/n he a fixed positive number with < A < r, 
let C, he the unique positive root of the equation 



and let 

Then, with k = \n an integer, one has, as n ^- oo, 



k 

log Fjn.riC) logC 

n 



In addition, a full expansion in descending powers of n exists. These esti- 
mates hold uniformly for A in any compact interval of [0, r]. 

Proof. One easily checks that the conditions of Theorem VIII. 8 in [8j are 
satisfied. D 

Fixing m and the fraction of errors k/mn = X/m, this allows us to control 
the asymptotics of pn^m,rik) for large n. 

Example 14. Figure [T] shows the probability Pn,m,r{k) that a word with k 
errors is correctible for Cm wr Sn for three different values of n, when m = 5 
(and thus r = [^^J =2). To reduce the error in the asymptotic formula 
for small values of k, we replace both the numerator En^m,rik) and denom- 
inator ("'™) oi pn,m,r{k) by the respective leading terms of the asymptotic 



expansion given by Proposition 13 (for the denominator we use r = m). One 
expects heuristically that first-order corrections to the leading asymptotics 
will largely cancel each other. Numerically this seems to be confirmed, as 
even for moderate values of n the agreement between the asymptotic result 
and the exact values is remarkably good. 

5.2. Asymptotics of Pn,m,r{p) for C^wr 5^ as tti — t- c«. Let us now 
consider the asymptotics of Pn,m,rip), given that the probability p of a single 
error is fixed. Here, we will discuss the case of m — t- oo as n is fixed. 

To deal with the truncated binomial sum Fm,r (x) , we will use the following 
integral formulation. 

Lemma 15. 

ds 



Fm,ri^) = Y^ f i'^ + SX) 



r+in 



where the contour of integration is a clockwise circle about the origin of 
radius less than one. 

Proof. Expand the integrand as a power series in s (which is absolutely 
convergent for |s| < 1) and integrate term-by-term. D 



10 
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0.2 0.3 

k/mn 



Figure 1. The probability Pn,m,r{k) that a word with k 
errors is correctible versus the error frequency k/mn for 
Cm^T^ Sn, when m = 5 (and thus r = [^y^J = 2). Shown is 
a comparison of the exact values (shown as points) and the 
asymptotic result (shown as curves) from Proposition 13 for 
n = 8, 16, and 32 (from right to left). 



As in the previous subsection, an asymptotic analysis of the integral in 



Lemma 15 is obtained from a saddle-point approximation. For m and r 



large, the behaviour of the integral is determined by the exponential of 

m log(l + sx) — r log s . 
There is a unique saddle a given by 



A. 
da 



log(l + ax) log a 

m 



namely, 



a 



X m — r 



Importantly, the saddle collides with the amplitude critical point s = 1 
when X = r/{m — r). This changes the asymptotic behaviour significantly, 
and we need a uniform asymptotic expansion to take this into considera- 
tion. The standard procedure here is to re-parameterise the contour by a 
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quadratic, i.e., 

logfl + sx) logs = -— - 7* + 5 , 

m I 

where 7 and b are determined by matching the location of the saddle point 
s = a with t = —7, and the location of the critical point s = 1 with t = 0. 

In the following result, erfc denotes the complementary error function, 
which is defined as 



We shall also need 

where 
and 



erfc(x) = —= / e ^ dt. 

p{P,x) = y/log{l + x) - /3\ogx - h{(3) 
h{l3) 
A{/3,x 



-/31og/3-(l-/3)log(l-/3), 
1 1 



J_ 



V2p{(3,x) 



Note that the radicand in p(/3, x) has a quadratic zero at x = /3/(l — /3), and 
that the correct sign has to be chosen to make p{/3, x) an analytic function 
near that point. 

Proposition 16. Let f3 = r/m he a fixed positive number with < /3 < 1. 
Then, with r = (3m an integer, we have, as m —)• 00, 



FmAx) = {I + XY 



^-eTfc{V^pi(3,x)) + ^Me-^P^'^'^^' 



(1 + 0(1)). 



In addition, a full expansion in descending powers of m exists. These esti- 
mates hold uniformly for /3 in any compact interval of [0, 1] and x in any 
compact domain of C\ M~ . 

Proof. This result follows from equation (9.4.22) in Section 9.4 of Bleistein 
and Handelsman [S] (with r = 0). □ 



As an alternative to using Lemma 15, one could have written Fm,r{x) in 
terms of an incomplete Beta function and used results of Temme [12] . 



The main result now follows immediately from Propositions 11 and 16 
Corollary 17. As m ^- c«, we have 
1 



Pn. 



\P) 



+ 



ei:ic{\/rnp{r/m,p/{l — p))) 
A{r/m,p/{l - p)) p(^/^^p/(i_p))2' 



V2 



mir 



(l + o(l)). 
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Figure 2. Shown are six curves: for ?n, = 8, 16 and 32 (from 
left to right), the respective curves for Pn,m,r{p) and for the 



asymptotic result from Proposition 17 are barely distinguish- 
able. 



Example 18. As we consider r = [ 
even for small values of m, such as m 



m—l 



2— J , asymptotically /3 = 1/2, but 
= 8 and /3 = 3/8, the expression in 



Proposition [16] provides a surprisingly accurate approximation, as can be 
seen from Figure [2} 



6. Conclusion 

In this note we have introduced a new algorithm for decoding the gen- 
eralised hyperoctahedral group Cmwr Sn- If n <C m, the performance of 
the algorithm is better from both the complexity perspective and from the 
error-correcting capability, when compared to the case n^ m. 

In particular, for large m the complexity of the algorithm is 0(mlog?7i), 
whereas it is 0{n^) 



11 



the 



for large n. As is evident from Proposition 
error-correcting capability is a monotonically decreasing function of n (this 
behaviour is demonstrated in Figure IT]). On the other hand, using the prop- 
erties of the complementary error function one can deduce from Corollary 
17 that for p < r/m the error-correcting capability increases monotonically 
as a function of m for m sufficiently large (this behaviour is demonstrated 
in Figure [2J) . 
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The results of the asymptotic analysis in Section [5] give reasonable approx- 
imations even for moderately small values of m and n. Extending Proposi- 



tion 16 it is possible to give refined asymptotic estimates which are uniform 
in m and n. 

As another direction, one could consider replacing the group Cm wr Sn 
with another wreath product, and modifying the algorithm to suit. First, 
one could replace the symmetric group Sn with another group K; however, 
this would give a much smaller number of codewords, and also the decoding 
algorithm would need to include a "membership-testing" algorithm (see Holt 
et al. [10]) to check whether the decoded permutation was an element of K. 
Second, one could replace the cyclic group with another group H; however, 
this would require a more sophisticated decoding process. 



Acknowledgements 

The authors would like to thank P. J. Cameron for numerous helpful com- 
ments and suggestions. The first author is a PIMS Postdoctoral Fellow at 
the University of Regina, and his research was supported in part by an EP- 
SRC CASE studentship, sponsored by the UK Government Communications 
Headquarters (GCHQ). 



References 

[1] R. F. Bailey, Uncoverings-by-bases for base-transitive permutation groups, Des. Codes 

Cryptogr. 41 (2006), 153-176. 
[2] R. F. Bailey, Error-correcting codes from permutation groups, Discrete Math, 309 

(2009), 4253-4265. 
[3] I. F. Blake, Permutation codes for discrete channels, IEEE Trans, Inform. Theory 

20 (1974), 138-140. 
[4] I. F. Blake, G. Cohen and M. Deza, Coding with permutations. Information and 

Control 43 (1979), 1-19. 
[5] N. Bleistein and R. A. Handelsman, Asymptotic Expansions of Integrals, Dover Pub- 
lications, Inc., New York, 1986. 
[6] P. J. Cameron, Permutation Croups, London Mathematical Society Student Texts 

(45), Cambridge University Press, Cambridge, 1999. 
[7] W. Chu, C. J. Colbourn and P. Dukes, Constructions for permutation codes in pow- 

erline communications, Des. Codes Cryptogr. 32 (2004), 51-64. 
[8] P. Flajolet and R. Sedgewick, Analytic Combinatorics, Cambridge University Press, 

Cambridge, 2009. 
[9] The GAP Group, CAP - Croups, Algorithms, and Programming, Version 4-4', 2004, 

|http: / /www .gap- system .org', 
[10] D. F. Holt, B. Eick and E. A. O'Brien, Handbook of Computational Group Theory, 

Chapman & Hall, Boca Raton, 2005. 
[11] S. Huczynska, Powerline communications and the 36 officers problem, Phil. Trans. 

Royal Soc. A 364 (2006), 3199-3214. 
[12] N. M. Temme, Uniform asymptotic expansions of the incomplete gamma functions 

and the incomplete beta function. Math. Comp. 29 (1975), 1109-1114. 



14 ROBERT F. BAILEY AND THOMAS PRELLBERG 

Department of Mathematics and Statistics, University of Regina, 3737 Was- 
CANA Parkway, Regina, Saskatchewan S4S 0A2, Canada 
E-mail address: robert.bailey@uregina.ca 

School of Mathematical Sciences, Queen Mary University of London, Mile 
End Road, London El 4NS, United Kingdom 
E-mail address: t.prellberg@qmul.ac.uk 



